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Abstract. Recent work of the operator algebraists P. Muhly and B. Solel, pri- 
marily motivated by the theory of operator algebras and mathematical physics, 
delineates a general abstract framework where system theory ideas appear in 
disguised form. These system-theory ingredients include: system matrix for 
an input /state/output linear system, Z-transform from a "time domain" to a 
"frequency domain" , and Z-transform of the output signal given by an obser- 
vation function applied to the initial condition plus a transfer function applied 
to the Z-transform of the input signal. Here we set down the definitions and 
main results for the general Muhly-Solel formalism and illustrate them for two 
specific types of multi-dimensional linear systems: (1) dissipative Fornasini- 
Marchesini state-space representations with transfer function equal to a holo- 
morphic operator- valued function on the unit ball in <C d , and (2) noncommuta- 
tive dissipative Fornasini-Marchsini linear systems with evolution along a free 
semigroup and with transfer function defined on the noncommutative ball of 
strict row contractions on a Hilbert space. 

1. Introduction 

In the last decade operator algebraists have worked on a sophisticated theory, 
based on the notions of C*- and ^"-correspondences, that encompasses many op- 
erator theoretical results, for instance, from characteristic model and interpolation 
theory [12[ lllj. both in time- varying and multivariable settings, as well as for a 
range of more exotic examples. In the recent paper |14j a transfer function realiza- 
tion result was obtained for the generalized Schur class (introduced in [13]) in this 
W* -correspondence setting; see also [4]. It is the purpose of this note to enrich the 
developed theory of [H| with a system theory interpretation, and show how the W*- 
correspondences can be specified to include some examples from multi-dimensional 
system theory. 

In Section [2] we delineate the two examples that we will obtain as special cases of 
the general Muhly-Solel framework: (1) commutative Fornasini-Marchesini systems 
with evolution along the positive-orthant integer lattice Z^, and (2) noncommuta- 
tive Fornasini-Marchesini systems with evolution along the free semigroup gener- 
ated by d letters. In Sections [3] and [4] we lay out the somewhat formidable array of 
definitions needed to formulate the general Muhly-Solel setup. The definition and 
basic properties of the abstract input /state/output linear systems arising in the 
general Muhly-Solel framework are given in Section[5J Section[5]is dedicated to the 
abstract point-evaluation and the associated abstract ^-transform, while in Section 
[7] we arrive at the frequency-domain version of the system and basic properties of 
the associated observation function and transfer function. An other type of lin- 
ear systems which fits the Muhly-Solel framework (1-D time-varying systems) and 
some other types which do not are mentioned in the final Section O Throughout 
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the exposition, new ideas and concepts are illustrated for the two basic examples 
as they are introduced. 

The notation is mostly standard, but we mention here a few conventions. For any 
index set Q we write £ 2 (fl) for Hilbert space of complex- valued square summable 
functions over f2: 

{/:n-C: E. e J/MI 2 <°°}, 
while denotes the space of bounded complex- valued functions over Q, For 

a Hilbert space H we use -^(fi) and £^{£1) as short-hands for £ 2 {Vt) ® 7i and 
£°°(n) (g) H, respectively. 



2. TWO EXAMPLES OF MULTIDIMENSIONAL LINEAR SYSTEMS AND THEIR 

TRANSFER FUNCTIONS 

In this section we review dissipative systems of Fornasini-Marchesini type; both 
commutative and noncommutative. In a somewhat more general form the noncom- 
mutative case was studied in [5j [6] with connections to robust control discussed in 
[7]. See also [3] for an overview of these types of systems and their applications. 

2.1. Commutative Fornasini-Marchesini systems. Consider a system evolv- 
ing over Z^ , the set of d-tuples of nonnegative integers, given by the system equa- 
tions: 



(2.1) 



x{n) = A\x{n — e{) + • • • + A d x(n — e d ) 

+Biu(n — ei) H h B d u(n - e d ), 

y(n) = Cx(n) + Du(n), 



where ei, . . . ,e d are the standard basis vectors in R d and we set x(n) = for all 
n £ Z d that are not in Z+, and with a system matrix of the form 



(2.2) 



Ai B 1 

A d B d 
C D 



' X ' 




' x d ' 


u 




. y 



Here Aj, Bj, C and D are bounded linear operators and X, U and y are Hilbert 
spaces. We further assume the system matrix (12.2[) to be contractive. In this case 
we speak of a dissipative system of commutative Fornasini-Marchasini type (see 
[9]). Note that the output sequence {y{n)} neZ ^ and state sequence {x(n)} neZ d^ 
are completely determined by the input sequences {u(n)} neZ d and the initial state 
x(Q). Moreover, the input, output and state sequences satisfy the energy balance 
inequalities: For each N G Z + we have 



E 



uj(n)\\y(n)\\ 2 

nGZ^;|)l|<JV 

< 



u(n)\\x(n)\\' 



• E 

neZl;|n|=iV 

^ uo(n) \\u{n) \\ 2 



< 



i|<jV 



where u(n) 
(nx, ...,n d ) 



— r with 



n>d and n\ — n\\ 



■ n d \ for each n = 



Zl. For a Hilbert space V we write £% M {^%) f° r tne weighted 
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space of sequences v = {v(n)} ne %d with v(n) <E U and 



H| 2 = J2 ^(«)ll«(r 



I 2 < oo. 



It then follows from the energy balance inequalities that if the input sequence u is in 
£u ,<* ( Z + ) ) the output sequence y is in l 2 y >w (Z^_ ) and 1 1 y \ | 2 < 1 1 u | | 2 + 1 1 x (0) 1 1 2 , while the 
state sequence x is bounded; more precisely sup ArgZ+ X)| n |=Af a -'( n )ll a; ( n )ll 2 < 00 ■ 
The state sequence is in t 2 x ^(Z^) if an additional stability assumption on state 
operators Ax, . . . ,A d is satisfied. 

We write M d for the open unit ball in C d : 

(2.3) M d := \z = {z l ,...,z d )&C l : ^|z,f<l 

and for a z = (z 1; . . . , z d ) in C d we use the notation z n for z™ 1 • • ■ 
n = (nx, . . . , rid) £ Given a bounded sequence v — {v(n)} neZ ^ , with v(n) in a 
Hilbert space V, we define a V- valued function v on B d via the Z-transform: 

(2.4) v{z) = ^ z n «(n). 

Provided that the input sequence is in the Z-transforms of the input, 

output and state sequences are related through the frequency domain equations: 

j x(z) = x(0) +(j2k=i z kA k ^jx(z) +(j2k=i z kB k ^u(z), 
{ y(z) = Cx(z) + Du(z), 
from which one finds that the input/inital state to output map is given by 

(2.5) y(z) = F{z)u{z) + W(z)x(0), 

where F and W are the transfer function and observation function given by 
F{z) = D + C(7-^Ax z d A d y x x 

(2.6) x(z 1 B 1 + --- + z rf 5 d ), 
W{z) = C(I-z 1 A 1 z d A d )-\ 

2.2. Noncommutative Fornasini-Marchesini systems. We write J- d for the 
free semigroup generated by the d-letter alphabet {1, ...,d}; elements of T d are 
words of the form a = ijv • • • *2*i with ij G {1, . . . , d} for j = 1, . . . , N. We 
also include the empty word that serves as the unit element with respect to the 
semigroup action: concatenation of words. With |a| we indicate the length of the 
word a. 

Now consider a system that evolves along the free semigroup T d given by the 
equations 

{x(k ■ a) = A k x{a) + Bku(a), k = 1, . . . , d, 
y(a) = Cx(a) + Du(a), 

where we have a system matrix of the form (|2.2[) that we assume to be contractive. 
For notational simplicity we assume U = y = C. Again, the output sequence 
{y(a)} a <zf d and state sequence {x(a)} a( zjr d are completely determined by the input 



(2.7) 
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sequence {u(a)} a eFd an d an initial condition: x(fi) = x®, and in this case we have 
energy balance inequalities of the form 



E b(«)ii 2 + E N«)H' 



< 



< E \\u(a)\\ 2 + \\xm\\ 2 

aeF d \\a\<N 

for each N 6 Z + . Hence an input sequence u from t 2 (!Fd) produces an output 
sequence y in l 2 (Td) so that ||?/|| 2 < ||u|| 2 + ||a;(0)|| 2 , and the state sequence x is 
bounded; more precisely sup ArgZ+ X)|q|=at 1 1 a? (a) 1 1 2 < oo. In case the state operators 
A\, . . . , Ad satisfy an additional stability the state sequence x is in i 2 x (Td)- 

For this example we shall use a frequency domain consisting of operator-tuples 
rather than tuples of complex numbers. We therefore assume that we are given 
an additional Hilbert space JC and a d-tuple T = (T±,...,Td) of operators on 
fC. Such a (i-tuple T is said to be a strict row contraction if the operator matrix 
[ Tj • • • Td 1 defines a strict contraction from /C d into /C. For a strict row 
contraction T we define the noncommutative functional calculus 

(2.8) T a =T iN ---T n if a = i N ---h, T % = I K . 

Then for W = C, X and id = {w(a)} a ^ d a bounded W-valued sequence we define 
a function w on the set of strict row contractions via the noncommutative Z- 
transform: For each strict row contraction T we set 

(2.9) w(T) = E w ( a ) ® Ta e W ® K )- 

Then the system (|2.7I) is described in the frequency domain by the equations 
( x(T) = x{%) + (J2 d k=1 A k ® T fc ) z(T)+ 

(2-io) +(Efc=i B *®r fc )«(T), 

1 y(T) = (C ® JjcJjbCT) + (D <8 Jjc)u(T), 
and the input/intial state-to-output map is given by 

(2.11) y(T) = F(T)u(T) + W{T)x{%), 

where we define the transfer function F and observation function W on the set of 
strict row contractions T by 

F(T) = D + C(I-^i®Ti Ad^Td)- 1 x 

x{B 1 ®T 1 + --- + B d ®T d ) 

d 

(2.12) = £> + E E GA a B k T a k 

k=l a^Fd 

W(T) = C{I-A X ®T X A d ®T d )- 1 

= ^{CA a )®T a . 

It is possible to symmetrize the noncommutative Fornasini-Marchesini system 
equations (|2.7[) to recover the commutative Fornasini-Marchesini system equations 
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(|2.1|) . To describe this in detail, we introduce the abelianization map 
(2.13) a(a) = [m, . . . ,n d ) (a E T d ) 

where, for j = 1, . . . , d, the integer nj is equal to the number of times the letter j ap- 
pears in the word a. Let {u(a), x(a), y(a) : a E Td} be a system trajectory satisfy- 
ing flMI), for each n E define u(n) = J2 a -. a(a)=n = a(a)=« 
F( n ) = S a; a (a)=ny( a )- Then one can check that the more detailed system equa- 
tions (|2.7p imply that the aggregate quantities {u(n),x(n),y(n) : n E satisfy 
the system equations (|2.ip . Moreover, in general if the ci-tuple T = (2\, . . . , T<j) of 
operators on JC is commutative, then T a = T a =: T n whenever a(a) = a(a') = 
n E Zl. and the value of the Z-transform w(T) in (|2.9|) can be expressed as 
w(T) = X)riGZ d w{n)T n . This applies in particular to the case where we take 
K, = C to be one-dimensional; then T = (T\, . . . ,Td) becomes a (i-tuple of scalars 
z = (zi, . . . , Zd) E B d , and the noncommutative Z-transform (|2.9[) collapses to the 
commutative Z-transform (|2.4[) (with W in place of v). 

3. C*- CORRESPONDENCE PRELIMINARIES 

Let A and £> be C7*-algebras and E a linear space. We say that E is an (.4, £>)- 
correspondence when E is a bi-module with respect to a right B-action and a left 
.A-action, and E is endowed with a B- valued inner product ( , )e satisfying the 
following axioms: For any A, \i E C, £, rj, ( E E, a G .4 and 6 G 6, it is the case that 

(1) (K + K,v)e = M^v)e + h(C,v)e\ 

(2) <£• b,r,} E = (Z,r,} E b; 

(3) (a • £,77) B = (£,a* • rj) E ; 

(4) «,t^ = 

(5) (£,t)E>0{mB); 

(6) (£, £) B = implies that f = 0, 

and that i? is a Banach space with respect to the norm 

IKUs = IK£.0*||J Ke£0. 

where || ||g denotes the norm of B. We also impose that 

(AO • b = £ • (Aft) and (Aa) • f = a • (Af ) 

for any A E C,a E A,b E B,£ E E. In practice we usually write ( , } and || 
for the inner product and norm on E, and in case A — B we say that E is an 
^-correspondence. 

Example 3.1. Part A: For the commutative Fornasini-Marchensini system inter- 
pretation, we take A — C{U) and E — colj =1 [£(U)] = C(U 1 U d ) 1 where U is a given 
Hilbert space. Then E can be seen as an .A-correspondence with left and right 
„4-action and inner product given by 

a ■ co\ d =1 [Tj] ■ a' = co\ d j=l [aT ja '], 
d 

(col^jT^col^Ji?,]} =^2^Tj. 

j=i 

Part B: For the noncommutative Fornasini-Marchensini system interpretation we 
take the same correspondence, but then with U = C. 
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Given two (.4, ^-correspondences E and F, the set of bounded linear operators 
from E to F is denoted by C(E,F). It may happen that a T G C(E,F) is not 
adjointable, i.e., it is not necessarily the case that 

<T£,7> = <£,r*7> (ZeE^EF) 
for some T* e £(F, E). We use C a (E, F) to indicate the set of operators in C(E, F) 
which are adjointable. As usual we have the abbreviations C(E) and C a {E) in case 
F = E. 

The third inner-product axiom implies that the left „4-action can be identified 
with a *-representation (meaning a nondegenerate *-homomorphism) p of A into 
the C*-algebra C a (E). In case this ^representation ip is specified we will occasion- 
ally write <p(a)t; instead of a ■ £. 

Furthermore, an operator T 6 £(E, F) is said to be a right B-module map if 

T(£ ■ 6) = T(£) ■ & (£eE,beB), 
and a left A-module map whenever 

T(a-£) = a-T(£) (fe£,o6>l). 

We will have a need for various constructions which create new correspondences 
out of given correspondences. 

Given two (A, /^-correspondences E and F, we define the direct-sum correspon- 
dence E(BF to be the direct-sum vector space E®F together with the left A- action 
and right ^-action and the direct-sum Z?-valued inner-product defined by setting 
for each £, f £ B, 7,7' £ F, a £ i and b G B: 

o-(€e7) = (o-0©(o-7), 

Ke7)-6 = K-&)e(7-6), 

(c © 7, c © 7')s®f = (e, ?}e + (7, y>F. 

Bounded linear operators between direct-sum correspondences admit operator ma- 
trix decompositions in precisely the same way as in the Hilbert space case (B = C) , 
while adjointability and the left and right module map property of such an operator 
corresponds to the operators in the decomposition being adjointable, or left or right 
module maps, respectively. 

Now suppose we are given three C* -algebras A, B and C, an (A, £>)-correspondcncc 
E and a (B, C)-correspondcncc F. Then we define the tensor-product correspondence 
E <S> F to be the completion of the linear span of all tensors £(87 (with £ G E and 
7 € F) subject to the identification 

(3.14) (e-6)®7 = C® (6-7) [ieE^ eF,beB), 

with left „4-action, right C-action and the C-valued inner-product defined, initially 

just for pure tensors, by setting for each G E, 7, 7' G F, a e ^4 and c G C: 

a • (£ <g> 7) = (a • £) <g> 7, (£ <g> 7) • c = £ <S> (7 • c), 
(£ ® 7, £' ® YW = <(£, Ob ' 7, 7')f, 
and then extending by linearity and continuity to E <g> F. In case the left action 
on F is given by the ^-representation p we occasionally emphasize this by writing 
E ® v F for E ® F '. It is more complicated to characterize the bounded linear 
operators between tensor-product correspondences. One way to construct such 
operators is as follows. Let E and E' be (A, ^-correspondences and F and F' 
(B, (^-correspondences. Furthermore, let X G C(E, E') be a right module map and 



A ^-CORRESPONDENCE APPROACH TO MULTI-DIMENSIONAL LINEAR DISSIPATIVE SYSTEMS 

Y € £(F, F r ) a left module map. We then write X ® Y for the tensor product 
operator in £(75 (g> F, E' <E> F') determined by 

(3.15) X <g> Y(f <g> 7) = ® (Y7) (C®76B®F). 

If, in addition, X is also a left module map, then X®Y is a left module map, while 
y also being a right module map guarantees that X ® Y is a right module map. 
Moreover, X ® Y is adjointable when X and Y are both adjointable, and in that 
case (X ® y)* = X* ® y* . Notice that the left action on 15 ® F can now be written 
asflu <p{p) ® If S C a (E ® F), where 7^ denotes the identity operator on F and 
(/? indicates the left „4-action on E. 

4. CORRESPONDENCE-REPRESENTATION PAIRS AND THEIR DUALS 

In the remainder of the paper we consider the situation where A = B, i.e, E 
is an ^-correspondence. We also restrict our attention to the case where A is a 
von Neumann algebra and E an _4-iy*-correspondence. This means that E is an 
^.-correspondence which is also self-dual in the sense that any right „4-module map 
p: E ^ Ais given by taking the inner-product against some element e p of E: 

(4.16) P (e) = (e,e p ) E eA. 

Suppose that we are also given an auxiliary Hilbert space £ and a ^-representation 
a: A — » £(£); as this will be the setting for much of the analysis to follow, we refer 
to such a pair (75, a) as a Correspondence-Representation pair, or CR-pair for short. 
We further assume that a is faithful (injective) and normal (er- weakly continuous). 
The Hilbert space £ becomes an (^4, (C)-correspondence with left ,4-action given by 
er: 

a ■ y — <j(a)y for all a 6 A and y 6 £. 

Thus we can form the tensor-product (A, (C)-correspondence E <g) CT £. We then 
write E a for the set of bounded linear operators /i : £ — > 75 (g>„. £ which are also left 
■4-module maps: 

(4.17) E a = {/j: £ -> E® a £: ficr(a) = (</j(a) <g> 7 £ ». 

It turns out that 75 CT is itself a iy*-correspondence, not over A but over the W*- 
algebra 

cr(^)' = {6 G £(£) : ba(a) = cr(a)6 for all a e .4} C £(£) 
(the commutant of the image er(_4) of the *-representation a in £(£)) with left and 
right er(_4)'-action and cr(A) '-valued inner-product (•, -)e° given by 

fi-b = fib, b ■ fi = (7e ®b)fi, v) E ° = v* [i. 

Let 1 : a (A)' — > £(£ ) be the identity ^-representation: (,(&) = & for each b S cr(^4)'. 
Then (i5 CT , i) is a CR-pair that is said to be dual to (75, cr); see Section 3 in [11] for 
further details. In particular, £ can be seen both as an [A, C)-correspondence (with 
left action given by a) and a (a(A)', C)-correspondence (with left action given by 

0- 

Example 4.1 (Continuation of Example 13. ip . 

Part A: We take £ = U and for the ^representation er the identity map: cr(T) = T 
for each T G C{U). Then E(g>£=U d , a (A)' = CI U ■= {^Iu : M e C} and 

75" = {co\ d ]=1 [^I u ] : lij £C}C C{U,U d )- 
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The spaces <j(A)' and E a shall be identified with C and C d , respectively. 

Part B: Let K, be some Hilbert space. We then take £ — K, and for a : C — > CQC) 

the ^-representation <r(X) = XI k. for each A S C. Then E (g>£ = JC d , a (A)' = £(£) 

and 

£ CT = {col^ljTj] : Tj E £(£)} = C(K,JC d ). 

In these two examples the duality mentioned above takes the following form. 
In Part A we start with E — C(U,U d ), A — C(U) and a the identity map on 
C(U). We then find the dual CR-pair {E a ,l) with E a = C d , a (A)' = C (after the 
identification) and t(X) = XIu- This is precisely the CR-pair we use in Part 2 (if we 
take K, = U). It then follows from Part 2 that forming the CR-pair dual to (E a , i) 
we return to (E, a). In this sense the two examples are also each others dual. 

Since E is an ^-correspondence, we can inductively define the self-tensor-product 
^-correspondences E® n via E® 1 — 1 and E® n+1 = E ® (E® n ) for n = 1, 2, . . . . 
Formally we set E®° — A. Using the notation in (|3. 15|) , the left ^-action on E® n 
can be written as <p(a) (g> I E ® n -i if we use ip : A — > C a (E) to indicate the left A- 
action of E. We may then also form the (A, (^-correspondences E® n ® a £. Besides 
the left .A-action, these Hilbert spaces also admit a left er(„4)'-action via the fit- 
representation b i — * I E ® n <g> b for b € <j(A)' . Furthermore, for a n € E a we can form 
the tensor product operator I E ®n ®/x that maps E® n ®£ into E® n+1 ®£ . In a similar 
way we define the a (A) '-correspondences < y E a )® n and (a (A)', (^-correspondences 
(E a )® n ®,£ for n E Z+. 

One of the consequences of the duality between (i?, cr) and (E a , t) is the following 
result from 

Proposition 4.2. For each n e Z + the map : {E a )® n ® u £ ~> E® n ® a £ induced 
by the identity 

®n(Pn <8> ■ ■ ■ ® Ml ® e ) = ® Mn) ' ' ' (Jfi ® M2)Ml e 

is a unitary left cr(,4)'-module map. 

Example 4.3 (Continuation of Example 14. ip . 

Part A: For n 6 Z + set d„ = #{a e ^d, |a| = n} = cP. Then 

I?®" =C(U,U dn ) and ® ff £ = U d ~ , 

while 

(^)®« = £(C,C rf ") = C d " and (E a )® n &,.£ =U dn . 

In particular, E® n <S> a £ and {E a )® n ® L £ are the same Hilbert space, and moreover, 
the left cr(A)' — C-actions on both spaces is just scalar multiplication. The left 
A = £(W)-action on U " is given by 

(4.18) T {®%e 3 ) = ®% 1 Te j , 

for T € C{U) and ej 6 U. Moreover, the map $„ defined in Proposition 14.21 is the 
identity map on U dn . 

Part B: Due to the duality between the two examples, as observed in Example 
14.11 we find: 

E ®n = c ^ C d n j = C d n and E ®n ^ g _ £<i„ ; 

(£ ff )®" = £(/C,/C d ") and (£ ff )® n ® t £ = /C d '\ 
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with left a{A)' = £(/C)-action on /C d ™ given by (|iT8| . with T E £(£), ej E /C. 
Moreover, is again the identity map on /C d ™ . 



5. Linear dissipative discrete-time systems 



" A 


B ' 




' X ' 




' E a ® X ' 


C 


D 




£ 




£ 



Let (E,a) be a CR-pair as in Section [U A (linear) system associated with the 
pair (E, a) is a set {A, B, C, D, tt} consisting of four left cr(^4)'-module maps: 



(5.19) 



and a ^-representation tt : c(A)' — ► C(X) that provides a left cr(^4)'-action on the 
Hilbert space X, called the state space, making X into a (a(A)' , C)-correspondence. 
In this paper we are primarily interested in the case that the system matrix (|5. 19|) 
is contractive, in which case we speak of a dissipative system. For simplicity we 
only consider here the case that the input space and output space are the same, i.e., 
both equal to £; various techniques exist to extend to the case with different input 
and output spaces, cf., [14l |4], 

Given a system {^4, B, C, D, tt}, for each n E Z + we define the operators: 



An = I(E°)®n ® A, 
C n = $„(J( B c)»n ®C), 



where is as defined in Proposition 
following linear system of equations: 



(5.20) 



x(n + 1) = A n x(n) + B n u(n), 
y{n) = C n x{n) + D n u(n), 



D n = $„(/(£„)»„ ® £>)$*, 

and associate with {A, B,C, D,ir} the 



(n E Z H 



where u(n),y(n) E E® n ® £, and x(n) E (E a )® n (g) X. Note that an input se- 
quence {w(n)}„ e z + and initial state a;(0) uniquely determine the output sequence 
{y(n)} neZ+ and state sequence {x(n)}„ eZ+ . 



Example 5.1 (Continuation of Example [ 
Part A: Let {A, B, C, D, tt} be a dissipative system associated with the CR-pair 
of Example 14.11 Part A. Then tt : C — > C(X) must be equal to 7r(A) = \Ix, and 
thus E v ® X — X d . The system then has the form 



(5.21) 



.4 
C 



B 
D 



' A 1 




A d 


B d 


C 


D 



' X ' 




' x d ' 


u 




u 



and, since cr(A) 1 = C, the CT(^l)'-module map properties are satisfied trivially. The 
input and output sequences u = {u{n)} n< zi + and y = {?/(") }nez + have elements 
u(n) and y(n) in U kn , while the elements of the state sequence x = {a;(n)} rie z + are 
in X kn . The system equations take the form 

{x(n + 1) = A n x(n) + B n u(n), 
y(n) = C n x(n) + D n u{n), 



(5.22) 



(n E 1 A 
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with 



diag* 



At ' 

A d 



B n 



diagfei 



/ 


" B 1 ' 




V 




) 



(5.24) 



diag^CC), D„ = diag t t 1 (I?) ) 

where diag forms a block diagonal operator, i.e., A n G C{X dn , X dn+1 ), B n 6 
Ci]A dn , X dn+1 ), etc. After making the identification of an element u(n) in U dn 
with a set of elements of U of the form {u(a) : a G J-d, \cx\ = n} via the bijective 
map 

(5.23) v n : {a G T d : \a\ = n) —> {1, . . . , d n }, 

v„(a) = 1 + X)™=i(i? _ 1)^ _1 if a = i n --'ii, 

and similarly for A'''" , and after unpacking the system equations (|5.22j) , we arrive 
at the noncommutative system 

x(k ■ a) = Akx(a) + Bku(a), k = 1, . . . , d, 
y(a) = Cx(a) + Du(a), 

that evolves over We obtain the commutative system (|2.ip after applying the 
abelianization procedure. 

Part B: In this case the system matrix of a dissipative system {A, B, C, D, n} 
associated with the CR-pair of Example 14. li Part B, is as in (|5.2ip . but with 
IA = /C, and tt is a ^-representation from £(JC) into C(X). Moreover, the operators 
in {A, B, C, D, it} satisfy the following intertwining relations: 

n(T)Aj = ^vr(T), n(T)B 3 = B,T, 

TC = Cn(T), TD = DT, 

for any T G C-{K,) and j — 1, .... d. In particular, D is a scalar multiple of 1%. 
Assume that tt is a scalar multiple of the identity representation; i.e., there exists 
a Hilbert space X so that X = X ® K, and tt(T) = <g> T for each T G £(/C). In 
fact, if K. is finite dimensional, then any representation 7r has this form (see [2]). 
Then (f5T2"5)) simplifies to 



(5.25) 



(5.26) 



As 



Aj ® / K , =Bj®I K , C = C® Ik, 



with A,- G G C(C,X) and C G £(A\C) for j = l,...,d. Due to this 

special form of the system operators, we see that, if the input signal {u(a)} Q g^ [i , 
identified with the actual input signal {u(n)}„ e z + as in Part A, and initial state x(0) 
are restricted to have the form {u(a)ko} ae j^ d and 3?(0)<8>fco, respectively, where u(a) 
is scalar- valued, x(0) G X and A:o is a fixed vector in JC, then the output sequence 
{y(ct)} ae jr d — {y{a)k^} ae jr d will have the same form as the input sequence, i.e., 
y{ot) G C, while the state sequence is of the form {x(a)} a ^jr d = {x(a) <8> fcojae^j 
with x(a) G X. From now on, we shall restrict the inputs and initial states to 
vectors of this form for a fixed vector ko G K,. In that case, the transformation to 
the noncommutative system equations occurs in the same way as in Part A, and 
we arrive at (|2.7p with Aj, Bj and C replaced with Aj, Bj and C. 

We now define F 2 (E) — ©^o-E®" - , which is an ^-correspondence with left A- 
action given by the ^representation ip ^ : A — > C a {T 2 (E)) that is defined via 

foc(a) = diag neZ+ ((fi(a) <g> J E ®„-i) 
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where if is used to indicate the ^-representation that provides the left ^-action on 
E. Then we can form the tensor product space T 2 (E,a) = T*{E) ® a £, which 
in turn is an (A, C)-correspondence. We refer to T 2 (E) and T 2 {E,<j) as the Fock 
spaces associated with the correspondence E and the CR-pair (E,o~), respectively. 

Besides the left .A-action, the Hilbert space !F 2 (E, a) is also equipped with a left 
<r(A)'-action that is given by the ^-representation b \— > Ij^/e) ® b £ C(J- 2 (E,a)) 
for b £ a(A)'. 

In a similar way as above, we define the Fock spaces T 2 (E a ) = (B^oiE^)®" 
and J r2 (E a ,Tr) = ^{E 17 ) ®^ X. Moreover, with T co {E' j ,tt) we denote the space 
of bounded sequences {x(n)}„ e z + with x(n) G ( K E (y )® n ® v X. 

Example 5.2. (Continuation of Examvle \5.1}) 

Part A: As observed in Example Ol E® n <& a £ = U dn , which, via the map v n in 
(|5.23|) . can be identified with ® a ef d ,\a\=rM- It thus follows that T 2 {E,a) can be 
identified with ^(Td), and similarly J 72 (E a ,n) with £ 2 x {T d ). For both spaces the 
left C-action is just scalar multiplication. Via the same identification, T°° (E a , it) 
can be seen as the subspace of fxi^d) consisting of elements {x(a)} a( zr d from 
2x(Fd) such that su PjVeZ+ J2\ a \=N \\ x ( a )\\ 2 < °°. 

Part B: Similarly as in part A, T 2 (E,a) and J r2 (E (7 ,Tr) can be identified with 
l 2 K (T d ) and^(^d), respectively, and T 00 (E a , it) with a subspace of £^{T d ). How- 
ever, in this case we have a left £(/C)-action that is given by entry-wise multiplica- 
tion with T G C(K) (resp. tt(T)) on the left. 

From here on we shall identify the spaces T 2 (E,o), T 2 {E a , tt) and T x '{E a ', tt) , in 
both parts A and B, with the appropriate £ 2 and £°° spaces over 

With the state operator A and two nonnegative integers n, m G Z + we associate 
the following generalized power. 

(A n ) m = (I {E „ )9m -i ® A n ) ■ ■ ■ {I E ° ® A n )A n , 

which defines an operator in C{(E a )® n ® £, [E a )® n+m <g> £). Here A° is to be 
understood as I[E<*)® n ®S- Note that (A n ) m — J( E ^j®n (g> A™ = (A" l )„. Thus we 
can unambiguously write A™ instead of (A n ) m or (A m )„. 

Theorem 5.3. Let y = {y{n)} n£ z + and x = {x(n)} n£ z + be the output and 
state sequences generated by the dissipative system (|5.20l) with input sequence 
u = {it(n)} nC z + and initial state x(0). Assume that u G J- 2 (E,a). Then y G 
T 2 {E,o) , x G T°°{E°,tt) and 

y = Tvu + Wvx(0), 

with T s G C{T 2 {E,a)) and W E G C(X, T 2 {E, a)) the left cr(„4)'-module maps 
given by 





Do 










' C A° - 




C^Bo 




••• 




CiAl 




C 2 A\B 


C 2 A° 2 B 1 


# 2 ... 




C2A 2 



Moreover, the operator matrix [ Ts Wy, ] is contractive, and thus ||y|| 2 < ||it|| 2 + 
||x(0)|| 2 . Finally, if, in addition, ||A|| < 1, then x G T 2 (E a ,tt). 

The fact that is a cr(^4)'-module map means that it commutes with the 
operators Ir 2 (E) ® b for each b G a(A)'. In addition, it is not hard to show (see 
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the proof of Theorem 3.5 in [14]) that Ts also commutes with I^(e) ® A* for each 
fi G E a . In order to make sense of I^{e) ® M as an operator on J- 2 (E, a), because 
H maps £ into E <S> £, we have to make the identification T 2 (E) ® E = T 2 [E\ 
These commutative relations are one way of defining the analytic Toeplitz algebra 
T+(E,<j) associated with the CR-pair (E,a) of |10J . There are various other ways 
to characterize 7+ (E, a) ; originally in [10] the space was generated by a collection 
of creation operators, while in [5] a Toeplitz structure characterization is given. 
The definition used here stems from [TT] . 

6. Point-evaluation 

Given a dissipative system {A, B,C, D,ir} as in (|5.19p . there is a notion of point 
evaluation giving rise to a ^-transform enabling the passage from the time domain 
to the frequency domain for this abstract setting. In the W-correspondence setting 
studied here, points are pairs (77, b) with 77 an element of the generalized disk 

D(E a *) := {n: 77* G E a , \\n\\ < 1}. 

and b G a{A)'. With a sequence / = {/(n)}„ eZ+ G T 2 (E, a) (so f(n) G E® n ® £), 
we associate a function / A on 0(E a *) x a (A)' with values in £ defined by 

00 

f A (v, b) = J2 V n (lE®r> ® b)f n ((77, 6) G B(E") x a(A)'). 

n=0 

Here we use rf 1 to indicate the n th generalized power of 77: 

??" = tj(Jb ® ??) • • • (Jb®»-i ® 77) : ® f -> £, 

where for n = we set 77° = Ig. Note that the generalized powers 77™ are well 
defined because 77 is a left A- module map; 77 inherits this property from the fact 
that 77* G E a is a left A- module map. The sum in the point evaluation is well 
defined because ||?7 ra || < ||?7|| n for each n. 

With a element x = {x(n)} n& z + G J-°°(E a ',ir) we also associate a function x A 
on H)(E' 7 *) x <r(A)' , but this point evaluation is essentially different from that in 
f 2 {E,o). Given /1 G .E CT , we define the operator : a (A)' — > -E" 7 by 

L^b = n®b = 11b (bea(A)'). 

Then L M is a right <7(„4)'-module map, and thus, for n = 1,2,3,..., we can define 
the n th generalized power: 

= {Lfj, O /( E <r)®n-l) • • • (L p ® I E °)E^ 

which is an operator from cr(A)' into (_E <T )®"; we set L° equal to the identity 
map on a(A)'. We further define L^^ — L M ® Ix and L™ j7r = L" ® Then 
one can compute that < IImII" f° r each n G Z + . In particular, for each 

?7 G D(i5 CT *) we have the operators L^* )7r and L™„ „.. Then, with a state sequence 
.x = {x(n)} n( zz + G J 700 (i?' 7 , 7r) we associate a function x A on IJ(iJ cr *) x ct(.A)' with 
values in X that is defined by 

00 

x A (rj,b) = ^(L; l , 7r )*(6«)/ (E . ) 0^i^)x(n), 

n=0 

for each (77, b) G ~B(E (T *) x cr(^4)'. In a different context, and restricted to the case 
x = {x(n)} ne z + G T 2 (E° , 7r), this form of point evaluation was first studied in [15j : 
see also [5]. 
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Example 6.1 (Continuation of Example 15. 2 j) . 

Part A: In this case B(i? cr *) can be identified with the unit ball B d defined in 
(|2.3[) (formally elements of D(£ ,<7 *) are row vectors) and a(A)' = C. Let r) — 
(zi, . . . , Zd) G M d . Then the generalized power rf 1 corresponds to the vector in C d ™ 
whose j entry is equal to 

z a = Zi n ■ ■ ■ z H if a := v n (j) =i n ---i 1 . 

Since z\,...,z& are scalars, and thus commute, the entries in rf 1 corresponding 
to ai and 012 (via the map v n ) that are mapped to the same element in lA_ via 
the abelianization map a in (|2.13[) are the same. Then for u — {u(a)} a £j^ d in 
T 2 (E : a) = £u(Fd), the point evaluation in a point (z, b) 6 M d x C is given by 

u(z, b) = z a bu{a) = b Y z n u{n) 

where u(n) — X) a (a)=n u ( a ) f° r eacn n e ^ n particular, the part of the point 
evaluation that comes from <j(A)' plays no real role, and will be left out in the 
sequel, i.e., we only consider b = 1 G C and shall write u(z) instead of u(z, 1). 
Using that T°° (E a , tt) can be identified with a subspace of l^(Td\ it follows in 
a similar way that the point-evaluation for an x = {x(a)} a fzr d in .T 72 (£? CT , 7r) is 
defined by: 

x(z) = Y z n Tc{n) (z e B d ), 
where x(n) = J2 a ( a )=n x ( a ) f° r eacn n e 

Part B: In the noncommutative setting D(i5 CT *) corresponds to the set of strict row 
contractions T = (Ti, . . . ,T d ) with Tj 6 £(/C) and cr(^)' = £(£). Similar to the 
case for Part A, the n th generalized power of T is defined via a noncommutative 
functional calculus: the j th entry of T n is given by T a in (|2.8|l in case a := v n {j)- 
Since the elements Tj of T do not (necessarily) commute pairwise, we obtain a 
noncommutative point evaluation for an element u = {u(a)} ae r d in F 2 (E,o) = 

u(T,b) = T a bu(a) 

ae^d 

where T is a strict row contraction and b 6 £(/C). The point evaluation of an 
element x = {x(a)} a£ r d in J 7 °°(E a ,n) = l 2 x {Td) goes through the ^-representation 
7r: 

*(T,6) = n(T a b)x(a) 

with again T a strict row contraction and b G £(/C). Now, under the restrictions 
of Example 15.11 we have input sequences of the form u = {u(a)ko} ae jr d , with 
u{a) G C, and state sequences of the form x = {x(a) ® fco} ae jr d , with S(a) G A", 
where fco is a fixed vector in 1C. In that case u(T, b) = ^2 ae:Fd u(a)T a bko and 
x(T,b) = J2ae^ d %( a ) ® (T a bko), and replacing fco with 6fco we see that we may 
restrict the point evaluation to the case that b = Ik.. 
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7. Frequency-domain equations 

With the point evaluations as defined in the previous section we have the follow- 
ing result. 

Theorem 7.1. Let {A, B, C, D, ir} be a dissipative system. Given an input se- 
quence u = {u(n)} n( zz + G J-^{E, a) and an initial state x(0), let y = {j/(n)}„ 6 z + G 
T 2 (E,a) and x = {x(n)}„ e z + G J 700 (E' T ,Tr) be the output and state sequences 
given by the system of equations (|5.20[) . Then the Z-transforms u A , x A and y A of 
u, x and y satisfy the frequency domain equations: 

x*( V ,b)=x{0)+L^Ax*( V ,b)+L^Bu*(ri,b), 

7.27) 

y A (r),b) = Cx A (r 1 ,b)+Du A (f 1 ,b), 

and we have 

(7.28) y A (r),b) =F b (»j)« a (» 7 ,6) + W e (»j)x(0), 

where the transfer function i 7 ^ and observation function W^, are given by 
F s (r ? ) = £» + C(7-L;, i7r 4)- 1 L;, i7r 5, 

7.29) 

w s (7 7 ) = c(/-L;, i7r A)- 1 . 



Example 7.2 (Continuation of Example 16. ip . 

Part A: For this example the system matrix has the form (|5.2ip . An element r\ of 
H)(E' 7 *) is identified with a point z = (z\, . . . , z d ) G B d , and in that case the operator 
L*.A G works out to be zi^i + • • • + z d A d . Similarly L*»B = Z\B\ + ■ ■ ■ z d B d 

and gives us 

F s (z) = D + C(I -z x A x ZiAj)-^ 

x (ziBx H h z d S d ) 

Ws(«) = C(/-«iA x z^)- 1 

for the transfer function and the observation function, in agreement with (|2.6[) , and 
fl7T2"g|) reduces to (|2~5)) . 

Part B: For this example, a(A)' — C(IC) is nontrivial, ir: u{A)' = C(K.) — > 
is a nontrivial representation, and D(i? '*) is identified with the set of strict row 
contractions T = [T x •■• T d ] G C{K. d ,K)- Nevertheless i?" 7 (gi^ <Y is again 
identified with «Y d and the system matrix has the form (15.2ip . Following the 
assumptions made in Example 15.11 i.e., X = X ® /C and the system operators 
are of the form (|5.26[) . if r] G H)(E cr *) is identified with the strict row contrac- 
tion T — [Ti • • • T^] , then it works out that the operators L**A G C{X) and 
L*,B G £(£, A?) can be identified with 

L*,A = A 1 ®T 1 + --- + A d ®T d , 

L*.B = B~!®Ti + ■ ■ ■ + B d ®T d 

whence the transfer function and observation function assume the form 

F S (T) = D + C ® I K {I -~Ai ®Ti AiSTdJ^x 

x (B 1 ®T 1 + ---B d (g>T d ), 



W S (T) = C ® 7 K (/ - Ax ® Tj A d ® T« 



'/J 



\-l 



in agreement with (|2.12p . and, moreover, (|7. 28[) collapses to (|2.1ip . 
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8. Other examples and other types of generalized systems 

Our main purpose is to illustrate the general Muhly-Solel framework of tensor 
algebras, generalized Schur class, and Hardy algebras arising from a CR pair for two 
specific cases, namely: (f) commutative Fornasini-Marchesini input/state/output 
systems and their associated transfer functions as holomorphic functions on the 
unit ball, and (2) noncommutative Fornasini-Marchesini systems with associated 
transfer functions defined on the noncommutative unit ball. We mention that 
the Muhly-Solel formalism has at least one other special case which is a familiar 
setting in system theory, specifically, time-varying systems with the Muhly-Solel 
Schur-class function identifiable as the input-output map for a discrete-time linear 
time- varying system. Details on this can be found in [4] . 

Still other types of systems and associated transfer functions have been in- 
troduced and studied recently. Specifically there is a system theory over hyper- 
holomorphic functions (with quaternions playing the role of the complex numbers) 
and a system theory for stochastic systems where the system parameters are them- 
selves random variables; in these instances the familiar algebra centering around 
pointwise multiplication of functions is replaced by something more exotic: Cauchy- 
Kovalevskaya product of hypcrholomorphic functions or Wick product of stochastic 
variables. An introductory survey to these matters can be found in pQ. 
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